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ABSTRACT 

A  new  procedure  is  described  for  determining  that  frequency  at  which  the  power  spectrum 
of  a  signal  has  its  absolute  peak.  The  salient  feature  of  the  procedure  is  that  it  does  not 
explicitly  involve  the  estimation  of  the  power  spectrum  of  the  signal  itself.  Specifically, 
it  is  shown  that  the  limit  of  the  iterated  normalized  autocorrelation  [see  Eqs.  (4)  and  (5)] 
of  a  function  f(t)  is  a  pure  cosine  wave  whose  frequency  corresponds  to  the  location  of  the 
maximum  energy  density  in  the  spectrum  of  f(t). 

Furthermore,  if  one  is  willing  to  accept  the  frequency  of  moximum  energy  to  within  a 
given  finite  spectral  resolution,  then  the  procedure  terminates  after  o  specified  finite 
number  of  iterations.  Results  from  a  computer  simulation  of  the  procedure  are  also  de¬ 
scribed. 

The  areas  of  application  of  this  procedure  are  discussed.  The  results  indicate  that  this 
method  of  detecting  a  signal  (i.e.,  by  the  peak  of  its  spectrum)  merits  further  considera¬ 
tion. 

It  is  important  to  note  that  the  effects  of  noise  have  not  been  considered  in  this  initial 
study;  the  results  apply  to  the  received  signol  only. 


This  technical  documentary  report  is  approved  for  distribution. 


l-ranklin  C.  Hudson,  Deputy  Chief 
Air  Force  Lincoln  Laboratory  Office 
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DETECTION  OF  THE  ENERGY  PEAK 
OF  AN  ARBITRARY  SIGNAL 


I.  INTRODUCTION 

In  communications  engineering,  it  is  often  useful  to  be  able  to  extract,  from  an  incoming 
signal,  that  frequency  which  contains  a  greater  energy  density  than  any  other  frequency.  In  this 
study,  a  new  procedure  is  described  for  determining  that  frequency;  the  significant  aspect  of  the 
procedure  is  that  the  power  spectrum  of  the  signal  need  not  be  calculated. 

The  fundamental  theorem  of  Sec.  Ill  describes  the  mathematically  interesting  result  in  the 
case  where  we  carry  our  procedure  to  the  limit.  The  more  useful  theorem  (Theorem  2,  Sec.  V), 
however,  describes  a  realizable  procedure  for  determining  the  frequency  of  maximum  energy 
to  within  a  finite  spectral  resolution. 

We  begin  by  defining  the  quantities  basic  to  the  procedure. 


n.  DEFINITIONS 

Consider  that  class  of  real  functions  or  signals  f(t)  whose  autocorrelation  function  g(t)  has 
the  following  properties: 

0  <  g(0)  <  «  (1) 

C  |g(t)  I  dt  <  »  (2) 

J  -OO 

g(t)  is  continuous  in  t  . 

Since  g(t)  is  an  autocorrelation,  it  is  therefore  an  even  function  of  its  argument.  We  now  define 
a  set  of  normalized  autocorrelation  functions  Rn(t)  as 


and 


R0(t)  = 


g(t) 

i(oy 


(4) 


Rn(t)  = 


Rn-l(T)Rn-l(t  +  T)dT 


/OO 

_o 


n  =  1,2,3,.. 


R  .(t)  dr 
n-1 


(5) 


We  recognize  that  R  (t)  is  the  normalized  autocorrelation  function  of  the  normalized  autocorre¬ 
lation  function  RQ(t),  etc.  We  may  thus  consider  Rn(t)  to  be  the  nth  iterated  normalized  auto¬ 
correlation  function  of  f(t). 
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Consider  the  Fourier  transform  Sn(a>),  defined  by 


SnM 


C  Rn(t)  dt 

—  GO 


(6) 


and  its  inverse 

Rn(t)  =  Sn(Cj)  eJ£l>t  ^  ‘  (?) 

Note  that  S  (01)  is  proportional  to  the  spectral  density  of  Rn_1(i);  in  particular,  Sq(oi)  is  propor¬ 
tional  to  the  spectral  densityt  of  f(t).  We  also  observe  that  Rn(t)  and  Sn(w)  are  both  real,  even 
functions  of  their  arguments. 

HI.  THE  FUNDAMENTAL  THEOREM 

In  Sec.  I  we  stated  that  our  procedure,  consisting  of  operations  limited  strictly  to  the  time 
domain,  enables  us  to  determine  the  peak  of  the  spectrum  of  a  function  f(t).  We  will  now  state 
the  theorem  involved. 

Theorem  1. 

If  S  (ca)  has  a  unique  absolute  maximum  at  0,  viz., 

SQ(0)  >  SQ(w)  for  a)  ^  ±©  ,  (8) 

then 

(9) 

lim  Rn(t)  =  cosQt 
n-*-°° 

and 

Sn(u)  |So(w)|2n  (10) 

sTJe)  =  lso(0)l 

Proof. 

We  first  demonstrate  that  exists  for  all  finite  n.  A  sufficient  condition1  for  the  ex¬ 

istence  of  Sn(co)  is 

C°°  ,  i  (11) 

\  |R(t)|dt<~  • 

V  -  oo 

We  will  establish  this  condition  by  an  inductive  proof  on  the  hypothesis  H,  which  states  that  R„(t) 

is  integrable  and  continuous  in  t.  Clearly,  by  Eqs.  (1)  through  (4),  we  have  Hq.  We  now  assume 

H  ,  and  show  that  this  implies  H  From  Eq.  (5)  we  may  form 
n-1  11 


r 


I  Rn(t)  I  dt  = 


si  I  /:. 


Rn-1(T)  Rn-1 


(t  +  t)  dT  I  dt 


(12) 


/: 


Rn-1(T)  dT 


t  Whether  this  is  the  energy  density  or  power  density  depends  upon  the  class  that  f(t)  falls  In.  See  Eqs.  (25) 
through  (28)  and  the  discussion  following  Eq.  (28). 

♦  See  Ref.  1,  p.  103. 
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But 


R 


At)  R  ,(t  +  t)  dr 
n-1  n-1 


Rn-t(T)|  |Rn-l(t+  T)|  dT 


thus 


R  (t) 

n  1 


dt  < 


[/r,  lRn-l(T)l  dTl 

/.:  Rn-1(T> 


(13) 


However,  by  assumption,  J"  |R  1(t)|  dT  <  °°.  Furthermore,  since  Rnl(t)  is  assumed  to  be 
continuous  and  Rnl(0)  =  1  (from  Eq.  (5)],  then  RZ_.,(t)  dr  >  0.  This  establishes  Eq.  (11). 
Now,  since  for  normalized  autocorrelation  functions^ 


R  „(t) 

n-1 


<  R  ,(0)  = 
n-1 


(14) 


we  may  write 

Rn-l(T)Rn-l(t+  T)^Rn-l<T> 


(15) 


Therefore,  we  may  bound  from  above  the  integrand  of  the  numerator  in  Eq.  (5)  by  a  function  which 
is  itself  integrable  and  continuous  (by  assumption  of  H  ^).  Thus,  Rn(t)  is  continuous*  This 
proves  H  and  completes  the  inductive  proof,  thus  establishing  the  existence  of  for  all 

fixed  n. 

We  now  use  the  fact  that  the  Fourier  transform  of  the  autocorrelation  of  a  real,  even  function 
is  equal  to  the  square  of  the  Fourier  transform  of  the  function  itself.  In  our  case,  we  have  al¬ 
ready  established  the  existence  of  the  Fourier  transform  of  the  function  Rn(t),  namely,  Sn(u>). 
Thus,  since  the  numerator  of  Eq.  (5)  represents  the  autocorrelation  of  the  function  Rn_^(t),  we 
transform  both  sides  of  Eq.  (5)  and  make  use  of  the  product  relation  above  to  obtain 


S  (w) 
n 


S  ,(u>) 

n-1 


r  R2  , (r)  dT 
n-1 


n  =  1,  2,  3,  .  .  . 


(16) 


We  now  concentrate  our  attention  on  the  frequency  9.  Forming  the  ratio  Sn(u))/Sn(9),  we  obtain 
from  Eq.  ( 16) 


S 

S 


n 


n 


(u>) 

(©) 


Sn-l(a,)lZ 

Sn-l(e)J 


It  is  then  clear  that 

Sn<“>  fSo<M>lZn 
Sn(0)  -  lSo(9)j 


(17) 


t  See  Ref.  1 ,  p.  61 . 
tSee  Ref.  2,  p.  67. 
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which  establishes  Eq.  (10).  We  develop  some  alternate  expressions  for  Sn(a>)  as  follows.  Let 


f®  R2  At) dr 

J_oo  n-1 


Then  Eq.  (16)  may  be  written  as 


where 


SnM  =  KnS;  (W) 


2  2 2  2n 
Kn  =  “n-1  an-2’  '  '  '  '  ao 


But,  by  Eqs.  (5)  and  (7), 


R  (0)  =  1  = 
n 


=  M" 

2*  I. 


S  ( co )  dco 
n 


Integrating  Eq.  (18)  and  applying  the  above  relation,  we  see  that 

K  -  *1 _  . 

n  „  ,n 

loo  So  <“>  d“ 

Thus,  we  obtain,  as  alternate  forms  for  Sn(w), 


S  ( w)  =  2 it 


2n 

So  <“> 

LI  s  f  on  a* 


S  (co)  =  27T 
ir 


tso(co)/s0(e)l 


11  ISo(0-)/So(e)]2  ^ 

We  now  proceed  to  establish  Eq.(9).  By  assumption 
S  (u>) 

°  Q.  <  1  for  all  uj 

so(0) 

Thus,  by  Eqs.  (10)  and  (20),  we  see  that 

Sn(  w)  j  0  cu  ±0 

11^(0)  =  |l  -=±0  . 

Furthermore,  it  is  clear  from  Eqs.  (7)  and  (14)  that  for  all  fixed  n. 


OO 

\  Sn(cu)  dcu  =  2irRn(0)  =  2ir 
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Since  the  integral  of  Sn(oj)  is  constant  and  since  Sn(cj)  is  vanishingly  small  compared  to  Sn(9)  at 
all  w  ^  ±6,  we  recognize  that  the  limit  at  n  -*  «  must  bet 


lim  S  (w)  =  7ruo(to  -  0)  +  ttuo(u;  +  6) 
n-**> 


The  transform  of  this  limit  function  (by  Eq.  7)  is 


lim  R  (t)  =  cos  0t 

n-*-oo 


which  establishes  Eq.  (9)  and  completes  the  proof  of  Theorem  i. 

Corollary. 

If  S  (u>)  has  a  finite  number  of  equal  absolute  maxima  at  the  frequencies  ©k  (k  =  1. 2,  .  .  .  ,K), 

then 

K 

lim  Rn(t)  =  -g  Yj  cos  ekl 


SnM  [So(w)l2 


Proof. 

The  proof  here  is  identical  with  that  of  Theorem  1  through  Eq.  (20)  (where  ©k  is  substituted 
for  O).  We  now  recognize  that  the  limit  of  Sn(u>)  must  be 


Sn<ek>  =  lSo(0k> 


lim  Sn(co)  =  £  Y  tuo(aJ  "  ek'  +  Uo(co  +  V1 
k=l 

The  transform  of  this  limit  is  clearly 

K 

lim  R  (t)  =  Y  cos  ®k^  ’ 


which  completes  the  proof  of  the  corollary. 

Let  us  now  consider  those  functions  f(t>  whose  autocorrelation  function  g(t)  has  the  proper  les 
listed  in  Eqs.(l)  to  (3).  We  define  functions  f(t)  to  be  of  finite  energy  if 


in  which  case 


0  <  C  |f(t)  |  dt  < 

J  —  oc 

GO 

g(t)  =  \  f(r)  f(t  +  r)  dT 

—  oo 


f  Where  u  (x)  is  a  unit  impulse  occurring 


ing  at  x  =  0. 
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Further,  we  define  f(t)  to  be  of  finite  average  power  if 


0  <  lim  ^  1^)1 2  dt  ° 
T-oo  L  1  J-T 


in  which  case 


g(t)  =  lim  Jp  C  f( 
T-»  J-T 


t)  f(t  +  t)  dT 


Functions  of  finite  energy  or  of  finite  average  power  therefore  satisfy  Eq.  (2),  and  if  their  auto¬ 
correlation  function  satisfies  Eqs.  (1)  and  (3),  our  results  hold*  We  note  here  that  for  signals  of 
finite  energy,  S  (o>)  represents  their  normalized  energy  density  spectrum,  whereas  for  signals 
of  finite  average  power,  Sq(u))  represents  their  normalized  power  density  spectrum. 

IV.  EXTENSION  TO  PERIODIC  AND  DISCRETE  TIME  FUNCTIONS 

In  order  to  extend  Theorem  1  (and  its  corollary)  to  periodic  and/or  discrete  time  functions, 
we  need  merely  to  replace  certain  integrals  in  Sec.  Ill  with  the  expressions  described  below. 

For  continuous  periodic  functions,  our  results  hold  if  we  redefine  the  limits  of  all  previous 
time  integrals  to  extend  over  a  single  period  (Tq,  for  example),  and  if  we  redefine  all  integrations 
with  respect  to  «  as  sums  over  the  discrete  set  of  harmonic  frequencies,  viz.. 


x(t)  dt 


.  rT  /2 

1  (  o' 


x(t)  dt 


y(«>  d« -*  X  y<wm} 


Zrrm 

“  rp 

m  1 


For  these  periodic  functions,  we  obtain  periodic  autocorrelation  functions  g(t)  where 

g(t)  =  TxT  CT°/2  f(T>  f(t  +  r)  dr  .  (31) 

ro  J-To/2 

In  addition,  from  Eqs.  (5),  (7),  and  (30), 

y  s  (W  >  =  1  (32) 

L  n  m 
m=-°° 

and  therefore  we  must  adjust  the  factor  of  proportionality  by  multiplying  the  right  side  of  Eqs.  (23) 
and  (24)  by  l/2rr. 


j-  The  properties  expressed  in  Eqs.  (0  and  (3)  are  most  eosily  stated  in  terms  of  g(t)  and  will  be  left  in  that  form. 
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For  discrete  aperiodic  time  functions,  our  results  hold  if  we  replace  all  time  integrals  with 
infinite  summations  over  the  discrete  time  variable,  and  also  redefine  the  limits  on  all  integra¬ 
tions  over  to  to  extend  over  the  finite  range  -(ir/At)  <  co  «  (ir/At),  where  At  is  the  uniform  incre¬ 
ment  between  adjacent  time  samples.  That  is. 


and 


OO 

x(t)  dt  -  Y  x(tm)  ' 
m=-°° 


pir/At 

y(co)  dco  —  At  V  y(co)  dco 

-ir/At 


(33) 


(34) 


whe  re 


t_  =  mAt 

m 


(35) 


For  these  discrete  functions,  we  obtain  discrete  autocorrelation  functions  g(tm).  Furthermore, 
we  find  that 


£ 


ir/At  2ir 

S(w)d«  =  -5f  . 

-?r/At 


(36) 


thus,  the  new  factor  of  proportionality  requires  that  we  multiply  the  right  hand  side  of  Eqs.  (23) 
and  (24)  by  l/At. 

For  discrete  functions  (of  increment  At)  which  are  periodic  with  period  PAt,  we  must  make 
the  following  changes:^ 

P-1 

r  x(t)  dt  -  y  x<v 

J  or,  .  „ 


(37) 


i=0 


P/2 


J=  \  y(w)  dco  -  ^  Y  y(«m) 

m=-P/2 


(38) 


whe  re 


t.  =  iAt 

l 


2irm 
“m  ~  PAt 

For  these  discrete  periodic  functions,  we  obtain  discrete  periodic  autocorrelation  functions  g(t±) 
Moreover,  we  have 


P/2 


l  S<"m>  =  P 
m=  -P/2 


(39) 


Therefore,  the  new  factor  of  proportionality  requires  that  we  multiply  the  right  hand  side  of 
Eqs.  (23)  and  (24)  by  P/2ir. 


t  For  convenience ,  we  assume  P  t©  be  an  even  integer.  When  P  is  odd  we  must  alter  the  limits  in  Eqs,  (38)  and 
(39)  slightly. 
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V.  FINITE  SPECTRAL  RESOLUTION— A  FEASIBLE  TIME  DOMAIN  PROCEDURE 

The  interesting  feature  of  our  iterated  aistocorrelaition  function  i*  that  one  may  determine 
the  energy  peak  (at  co  =  0,  for  example)  of  a  signal  with  operations  atrictly  in  the  time  domain. 

In  practice,  one  usually  has  a  signal  f<t)  of  finite  duration  (T)  to  consider  t  As  a  result,  the 
autocorrelation  of  the  signal  RQ(t)  will  be  zero  outside  the  interval  |t|  *  T.  However,  the  for¬ 
mal  procedure  for  calculating  Rn(t)  as  described  byEq.  (5)  indicates  thatRn<t)  will,  in  general,  be 
nonzero  in  the  interval  |t|  <  2nT.  From  a  practical  point  of  view,  this  requires  an  exponentially 
increasing  complexity  (in  either  equipment  or  computation).  At  the  same  time,  we  have  an  expo¬ 
nential  rate  of  convergence  to  our  limit  function  where  the  exponent  is  2n  [as  may  be  seen  from 
Eq.  (10)).  If  I*  clear  that  this  rate  of  convergence,  although  exponential,  is  nevertheless  depend¬ 
ent  upon  the  shape  of  SQ(w).  Note,  however,  that  the  resolution  is  (theoretically)  perfect,  i.e.,  we 

are  guaranteed  to  converge  exactly  on  the  value  0. 

The  price  for  perfect  resolution  is,  as  always,  extremely  high  and  one  which  we  are  not 
willing  (or  able)  to  pay;  namely,  that  we  require  an  unbounded  number  of  calculations  if  we 
insist  on  passing  to  the  limit  n  -  «  .  Two  separate  aspects  of  the  complexity  of  the  process 
grow  without  bound:  the  number  of  iterations  n,  and  the  range  |t|  <  2nT  in  which  Rn(t)  must  be 
calculated.  If  we  are  willing  to  sacrifice  some  resolution,  then  we  may  control  both  of  these 
quantities  as  follows. 

Let  f(t)  be  the  signal  whose  energy  peak  (at  u  =  ©)  we  desire,  where  we  assume  f(t)  to  be  of 
finite  duration  T  sec,  i.e.,  f(t)  =  0  for  t  <  0  and  t  >  T.  For  this  function,  we  make  the  additional 
assumptions  stated  in  Eqs.  (t)  through  (3).  Define  Rn(t)  and  Sn(w>  as  in  Eqs.  (4)  through  (7).  We 
introduce  the  periodic  function  f  (t)  of  period  Tq  ^  2T  such  that 

fp(t  +  kTQ)  =  f(t)  ,  0«t«To  ,  k  = - -2, -1,0,  1,2 .  (40) 

For  this  function,  we  define  sj^o^)  and  R^t).  as  in  Eqs.  (4)  through  (7),  with  the  changes 
indicated  in  Eqs.  (29)  through  (32);  the  superscript  P,  serves  to  distinguish  these  from  Sn(w) 
and  R  (t)  which  refer  to  f(t). 

Below,  in  Theorem  2,  we  show  that  an  estimate  §  of  the  energy  peak  in  fp(t)  is  a  close  ap¬ 
proximation  to  O  when  Tq  is  sufficiently  large.  In  addition,  we  describe  a  procedure  for  ob¬ 
taining  this  estimate  to  within  a  given  accuracy  after  a  finite  number  (nQ)  of  iterations.  Further¬ 
more,  since  periodic  functions  have  periodic  autocorrelation  functions,  we  recognize  that  Rr  (t) 
is  periodic  (with  period  T  )  for  all  n.  Consequently,  we  need  carry  out  our  calculations  of 
R(P)(t)  over  only  one  period.  Thus,  by  agreeing  to  estimate  0  within  a  slight  (arbitrarily  small) 
error,  we  have  been  able  to  reduce  our  system  to  one  which  goes  through  a  finite  number  gfjjL 
e rations  over  a  finite  interval  (0,  Tq)  which  is  fixed  with  respect  to  the  number  of  iterations.  All 
this,  of  course,  comes  about  by  allowing  an  uncertainty  of  size  (p,  let  us  say)  in  the  estimate  of 
the  frequency  of  maximum  energy  density. 

The  relationship  between  the  time  limited  signal  f(t)  and  its  periodic  counterpart  fp(t)  needs 
further  elaboration.  Indeed,  we  recognize  that  the  autocorrelation  functions  of  these  two  signals 
obey  an  inverse  of  the  sampling  theorem.  Specifically,  sinceTo^2r,  Rq  (t)  contains  all  the 
information  about  RQ(t);  consequently,  S(oP)(a,m)  must  be  proportional  to  SQ(w)  at  co  =  «m  =  27rm/T 
We  may  calculate  this  factor  of  proportionality  from  Eqs.  (6)  and  (29),  viz., 

t  Therefore,  the  signal  wilt  be  of  finite  energy  [see  Eq.  (25)]. 
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S(P)(u>  ) 

q  '  m 

So(wm>  ' 


/  T°/2  R<P>(t)  e 

T°  -V2 


-iu>  t 
J  m 


°  m  f  R  (t)  e  dt 

J-oo  o' 

The  integrand  in  the  numerator  disappears  in  the  interval  T  <  |t|  <  Tq/2  and  the  integrand  in  the 
denominator  disappears  over  the  interval  |t|  >  T.  Therefore,  since  (Tq/2)  ^  T,  we  find  that 


S!,P)<-m> 


T^So<"m> 

o 


Recognizing  that  these  two  spectra  are  proportional  is  essential  to  one’s  understanding  of 
Theorem  2.  A  typical  case  is  shown  in  Fig.  i,  where  the  continuous  spectrum  SQ(u>)  is  shown 
as  the  continuous  envelope,  and  the  line  spectrum  TQSj)P)(cum)  is  shown  as  vertical  lines. 


Fig.  1.  Comparison  bet-ween  S  (u)  and  S  (u  ). 


i 

y 

% 

y 

j 

LjJk.i  *  •  — 

||V 

' 

‘•v 

■•V 

t  j 

pw 

0  „ 

'  “m'2Tm/To 

Observe  that  the  distance  between  "samples"  of  Sq(lj)  is  27r/TQ.  We  define  an  integer  M  such 
that  and  u>M+1  surround  the  spectral  peak  at  w  =  0  (the  same  clearly  holds  true  for  -  wM> 

-0  and  -wM+1).  Further,  we  define  p  to  be  that  frequency  at  which  Sq(w)  has  its  second 

greatest  local  maximum. 

In  estimating  e,  two  separate  problems  confront  us:  first,  how  do  we  eliminate  competing 
peaks  (such  as  at  u  =  ±p);  and  second,  how  do  we  converge  to  the  neighborhood  of  6  within  the 
absolute  peak  itself.  The  first  problem  may  be  solved  by  insisting  that  the  spacing  be 

fine  enough  to  insure  that  at  least  one  sample  in  the  neighborhood  of  O  (either  cuM  or  cuM+1>  is 
greater  than  the  maximum  possible  sample  near  u  =  p.  Thus,  an  important  parameter  of  the 
spectrum  S  (w)  may  be  expressed  as  a  lower  bound  A,  for  the  difference  between  these  two  peaks, 

viz. , 

so(e)  -  so(p)  >  a  .  (42) 

In  considering  the  problem  of  convergence  in  the  neighborhood  of  0,  we  recognize  that  the  shape 
of  the  spectrum  in  this  region  is  crucial.  Specifically,  we  require  that  the  spacing  2tt/To  be 
small  enough  to  guarantee  that  some  lines  in  the  discrete  spectrum  fall  within  the  range  of  the 
dominant  peak  (the  narrower  the  peak,  the  finer  must  be  the  spacing).  On  the  other  hand,  too 
many  spectral  lines  in  the  immediate  vicinity  of  9  will  result  in  a  large  number  of  required 
iterations,  because  the  largest  sample  will  not  differ  significantly  from  its  neighbors.  We  may 
discuss  the  width  or  sharpness  of  the  spectrum  in  the  vicinity  of  its  absolute  peak  by  considering 
the  curvature  of  Sq(w)  in  terms  of  the  magnitude  of  its  second  derivative  d  SQ(w)/dw  .  When  the 
second  derivative *is  large  (in  magnitude),  then  the  peak  is  narrow,  and  vice  versa.  Thus,  we 
are  led  to  consider  upper  and  lower  bounds  for  this  quantity  in  the  neighborhood  of  0,  viz.. 
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(43) 


d  S  (a;)  ,  47T 

-A<  - 2^—  <-a  for  |o>  -  0|  <  Tp— 

du>  0 

The  preceding  discussion  deals  qualitatively  with  those  factors  that  determine  the  required 
number  of  iterations  and  the  spacing  of  samples  in  the  discrete  spectrum.^  In  addition  to  these 
considerations,  the  sampling  must  be  fine  enough  to  satisfy  the  spectral  resolution  requirements 
of  the  user  of  these  results.  Stated  precisely,  we  have  the  following  theorem. 

Theorem  2. 

Consider  any  function  f(t)  of  duration  T  sec,  with  a  unique  absolute  energy  peak  at  w  =  ±0, 
whose  energy  density  spectrum  S o(o>)  satisfies  the  conditions  of  Eqs.  (42)  and  (43)  with  respect 
to  the  three  positive  parameters  (a.  A,  A).  Then,  for  a  given  required  frequency  resolution  p, 
there  exists  a  procedure  (defined  below)  which  will  calculate  a  number  0  after  nQ  iterations  of 
the  autocorrelation  of  fp(t)  (see  Eq.  (40)]  over  the  range  |t|  ^  Tq  such  that 

I  0  -  0  I  £  P 


where 


n 


o 


l°g2 


(44) 


and 


€ 

O 


a 


and 


T  =  max  (tt 
o 


'A  t-  4a 

2a 


2T ) 


(45) 


(46) 


✓N  {  P) 

We  begin  by  describing  the  procedure  by  which  we  determine  the  estimate  O  from  Rn  (t). 
We  recognize  that  if  one  line  (say  at  -  0  )  of  the  spectrum  of  the  periodic  function 
R(P)(t)  is  "sufficiently"  large*  compared  to  all  the  other  lines,  then  Rr  (t)  will  appear  as  a 
"noisy"  cosine  wave  at  frequency  0  rad/sec.  Our  procedure,  then,  is  to  count  the  number  of 
"noisy"  zero-crossings  of  this  function  for  a  known  time  interval;  if  we  are  successful  in 
counting  only  the  true  zero-crossings  of  the  pure  cosine  wave,  we  can  then  ascertain  0p  exactly. 
We  define  a  noisy  zero-crossing  counter  as  follows.  Consider  a  6 -threshold  detector  with  hys¬ 
teresis  defined  by  the  transfer  characteristic  shown  in  Fig.  2.  Further,  define: 

Z(t)  =  number  of  zeros  (or  counts)  recorded  by  the 
noisy  zero-crossing  counter  in  t  sec 


and 


Z(0)=  1 


tThe  sample  spacing  may  be  adjusted  by  changing  the  value  of  Tq 


u  —  u  . ! 
I  m  m-t 


=  2it/T 


tit  is  clear  that  by  increasing  T0,  0  can  be  made  as  close  to  0  as  one  desires;  also,  by  increasing  the  numbet 
of  iterations,  s(P>(6  )  can  be  made  Arbitrarily  large  compared  to  all  other  spectral  components  (see  Theorem  1). 
n  p 
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The  noisy  zero -crossing  counter  consists  of  a  8 -threshold  detector  with  hysteresis  followed  by 
a  simple  counter  which  registers  a  count  each  time  the  detector  changes  state  (in  either  direction) 
as  shown  in  Fig.  3.  It  will  be  shown  later  that  a  setting  of  6  =  i/9  defines  a  noisy  zero-crossing 
counter  which  produces  an  estimate  0  consistent  with  Theorem  2,  that  is,  <5n  ^  6  =  1/9  for  n  Ss  nQ 
[where  6n  is  defined  by  (*S)  below]. 


♦output 


If  i«eur 


8-threshold 

OETECTOR  WITH 
HYSTERESIS 


STATE 

CHANGE 

COUNTER 


Fig.  2.  Transfer  characteristic  of  tha 
5-threshold  detector  with  hysteresis. 


Fig.  3.  The  noisy  zero-crossing  counter. 


With  no  loss  of  generality,  we  may  write 

R^P)(t)  =  Bn  cos  Opt  +  bn(t)  <47> 

This  form  exposes  the  pure  cosine  wave  that  we  wish  to  detect,  and  groups  the  transform  of  all 
the  other  spectral  components  into  the  function  bn<t).  Define  6n  by 


Then,  if 


lbn(t)l<6n  for  1*1  <T0 


B  i  26 
n  n 


we  are  guaranteed  that  Z(t)  will  increase  by  unity  each  time  Bn  cos  ept  passes  through  zero, 
as  may  be  seen  in  Fig.  4.  That  is,  when  Bn  cos  0pt  is  negative  R^tt)  <  and  when  Bn 
cos©  t  is  positive  R<P>(t)  >-  6n,  thus  insuring  that  a  single  zero-crossing  (which  is  equivalent 
to  a  s'ign  change)  of  B^  cos  ©  t  cannot  generate  more  than  one  count.  Furthermore,  by  insisting 
that  Bn  >  26  ,  we  guarantee  that  each  time  Bn  cos  Opt  passes  through  zero,  we  must  get  at 
least  one  coSnt.  Thus,  Z(t)  -  1  will  count  the  number  of  zero-crossings  of  Bn  cos  0pt.  Now, 
since  Tq  is  a  multiple  of  the  period  of  Bn  cos  0pt,  we  determine  0p  from 

e  . 

P  Tcs 


b„  cos  U-n-iiTfl 


Fig.  4.  Description  of  the  reiotion  between  B^  and  8^ 


THIS  DISTANCE  MUST 
BE  AT  LEAST  8„ 
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Generally,  if  we  observe  Z(t)  for  a  time  r.  where  0  <  t  <  Tq,  then  we  may  ask  how  large  an 
error  is  made  in  our  estimation  of  6p;  in  particular,  how  large  is  |  6p  -  [Z(t)/t]  it  \ .  To  answer 
this,  we  use  Euclid's  algorithm  to  express 


t  ,  _r_ 

P  -  q  +  p 

where  P  is  the  period  of  cos  0  t  (i.e.,  P  =  2r/ep)  and  q  and  r  are,  respectively,  the  quotient 
and  remainder  of  r/P,  where,  of  course,  r  <  P.  But,  q  =  [Z(r)  -  1]  A  since  we  increment  Z(t) 
twice  for  each  period  P  of  cos  St.  Thus, 

r  _  Z(r)  -  1  _r_ 

P-  "  2  P 

or,  multiplying  by  27r/r,  we  obtain 

2 ir  _  ttZ(t)  x  7 r  „  £  , 

-p- - t~  +  T  P  ■ 

Now,  since  0p  =  2tt/P,  and  since  r  <  P,  we  have 

I e  ■  (5il 

1  p  T  T 

/p\ 

Thus  we  conclude  that  the  noisy  zero-crossing  counter  determines  from  Rn  (t)  the  frequency 

O  to  within  the  accuracy  described  by  Eq.  (51).  The  error  in  this  determination  may  be  made 

P 

arbitrarily  small  by  increasing  t. 

We  now  obtain  bounds  for  the  magnitude  of  the  difference  between  spectral  components  in 
the  neighborhood  of  0.  By  Taylor's  development,  we  have 

dS0(-’)  t  ,.,-co’)2  dZSQ(g) 

S0(")=  S0(co’)+  (co-  co')  — aiT—  +  —  du)2  ’ 

where  «r  lies  between  co  and  a;1.  Furthermore,  expanding  dSo(co')/dco  we  have 


dS  (to')  dS  (0)  d  S  (a) 

— ° -  =  -2 -  +  (co-  -  0)  - V" 

dco  do;  dco 


where  a  lies  between  co’  and  0.  Combining  these  two  equations,  and  recognizing  that  the  slope 

of  S  (co)  is  zero  for  co  =  ©,  we  obtain 
o 

d  So(a)  f co  -  co')Z  d  SQ(J)  (52) 

So(co)-So(co')  =  (co-co-)(co'-G)  — 

For  future  reference,  we  now  establish  a  lower  bound  for  this  difference  in  the  case  co '  =  «M* 
co  =  coM_r  Recall  that  -  «M_4  =  2*/T0  where  coM  is  defined  as  shown  in  Fig.  1 .  Thus, 

Eq.  (52)  becomes 


Q)  dw2 


d2So(u)  2xZ  d2So(u) 


Since  <r  and  a  both  lie  within  the  range  0  ±  47r/TQ,  we  may  apply  Eq.  (43)  and  in  obtaining  a 
lower  bound  for  the  above  expression,  we  may  set  =  ©,  yielding 
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(53) 


lSo("M)“So(wM-l)l  > 


Note  that  a  similar  equation  holds  for  o>M+1  and  “^+2' 

We  now  consider  the  conditions  necessary  to  insure  that  the  sample  at  either  a>M  or 
exceeds  the  maximum  possible  sample  at  u>  =  p.  Recognizing  that  one  of  the  two  samples  sur¬ 
rounding  6  must  lie  within  */T  of  0,  we  apply  Eq.  (52)  [with  «’  =  6,  a;  =  0  *  (ir/TQ)]  and  obtain 


Sq(0  *Tjr-)-So(0)|  = 


t2  d  So(u) 


Application  of  Eq.  (43)  gives  us 


|S  (e:k^)_so(e)U-iX  . 

o  £  l 

o 

(P) 

Thus,  the  minimum  value  of  the  maximum  sample  of  Sq  (<^m>  near  0  is 

where  we  have  made  use  of  Eq.  (41).  We  choose  to  require,  at  this  point,  that  this  minimum 
value  should  exceed  SQ(p)/To  by  at  least  the  quantity1-  e  2a  (tr  /Tq).  Thus, 

Applying  Eq.  (42),  we  obtain,  as  our  first  condition  on  Tq  [see  Eq.  (46)], 


To»vir  ■ 

We  have  eliminated  the  problem  of  competing  peaks. 

We  may  now  discuss  convergence  of  Rj[PNt),  in  the  neighborhood  of  9.  theorem  1  states* 

that  the  iterated  normalized  autocorrelation  (t)  of  fp(t)  converges  to  cos  0pt  where  SQ  (ep)> 

)  for  all  u  4*0.  However,  as  stated  in  the  corollary  to  Theorem  1,  if  there  are  a 
o  m  m  p  /pi 

number  of  equal  absolute  maxima,  then  R^  '(t)  will  converge  to  a  sum  of  cosines  at  the  various 
frequencies  of  the  equal  peaks.  We  have  postulated  that  So(u>)  has  a  unique  peak  at  u  -  O.  Never¬ 
theless,  SoPNu>m)  may  have  two  equal  peaks  at  =  <uM,  “M  +  1  as  shown  in  Fig.  5.  On  the  other 
hand,  S*P*(uM)  and  S^P1(w^+1)  may  be  arbitrarily  close  in  magnitude,  depending  upon  SQ(ti>)  and 
the  frequency  sampling  rate  2x/Tq.  In  order  to  handle  this  annoying  circumstance,  we  express 
R^P*(t)  (again,  with  no  loss  of  generality)  as 

R<P)(t)  =  C  1/3  cos  M  M  +  (1  -  fi)  cos  (M  +  1)1  +  cn(t)  .  (56) 

"  nl  Xo  o  i 


tNote  thot  t  is  just  the  lower  bound  in  Eq.  (53). 
o 

$  We  use  this  theorem  with  its  extension  to  periodic  functions;  see  Sec.  IV. 
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(  P) 

This  form  for  R^P\t)  exposes  the  interaction  of  the  two  largest  components  of  anci  groups 

the  transform  of  all  the  other  components  into  a  single  function  cn(t).  The  variable  0  allows  us 


Fig.  5.  Two  equal  peaks  in  S 


(P), 


“  )• 
m 


to  investigate  the  effect  of  the  relative  magnitudes  of  these  two  components,  where,  of  course 
0  4P  4  1.  It  is  clear  that  cn(t)  is  a  sum  of  cosines,  each  with  zero  phase;1'  therefore. 


cn(0)  *  ,Cn(t)l  ' 

Furthermore,  since  R^P\o)  =  1,  we  note  from  Eq.  (56)  that 

C  =  1  -  c  (0) 
n  n 


(57) 


(58) 


We  concentrate  our  attention  on  the  function 

e(/3,  t)  =  0  cos  ^  M  +  <1  -0)  cos  ^  (M  +  1)  (59) 

o  o 

and  investigate  the  behavior  of  its  maxima  and  minima  in  order  to  determine  the  effect  of  passing 
R(P)(t),  as  expressed  in  Eq.  (56),  through  the  noisy  zero-crossing  counter.^ 
n  Define  SR  to  be  that  instant  when  the  time  derivative  of  e(0,t)  has  its  k  zero,  i.e., 


de(0,S  ) 

“TT*--0  . 

Lemma  1* 

We  make  use  of  Lemma  1  as  follows.  For  all  0  <  0  .£  1 ,  and  for  all  k  and  t  such  that 

0«Sk^Uk.<T<To/2  , 


then 


|e({  ,r)|  <  |e((J,Sk)|  . 

Application  of  this  lemma  will  be  made  when  we  consider  Eq.  (56),  namely. 
R^P)(t)  =  Cne(/3,t)  +  cn(t) 

and  compare  it  to  Eq.  (47), 

R^P)(t)  =  Bn  cos  Opt  +  bn(t)  . 


(60) 


t Recall  that  S  (u)  is  always  real  and  even, 
n 

t  See  the  Appendix  far  proof  of  lemmo  1 . 
i  See  the  proof  of  Lemma  1  for  a  definition  of 
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Whereas  e(0,t)  is  not  a  pure  cosine,  it  does  have  the  same  number  of  relative  peaks  (maxima 
and  minima)  as  some  cosine  of  frequency  (<u',  for  example)  where 


uM4  u>'  4 


(61) 


The  same  can  be  said  about  the  number  of  zeros  of  e(£,t).  We  may  now  write  an  equation,  similar 
to  Eq,  (49),  which  places  a  condition  on  the  magnitude  of  the  relative  peaks  of  e(/3,t)  over  a  range 
of  t.  Specifically,  if  we  expect  the  noisy  zero-crossing  counter  to  calculate  the  frequency  a)' 
when  R^PVt)  (as  expressed  in  Eq.  (56)1  is  present,  then  we  require  [see  the  proof  of  Eq.  (49)]  , 

that 

Cn  |e<0,Sk)|  >2  |cn(t)|  . 


By  Lemma  1,  therefore,  we  require  that 

cn  |e(f  ,r)|  ^2  |cn<t)| 

for  0  <  Sk  ^  uk  ^  t  <  Tq/2.  By  Eqs.  (57)  and  (58)  we  obtain,  as  an  equivalent  condition. 


|e(|  ,  t)  |  >2 


cn(0) 

1  -  cn(0) 


(62) 


Thus,  we  propose  to  choose  some  r  in  the  range!  0  <  t  <  (Tq/2)  and  agree  to  observe  Z(t)  only 
Up  to  t  =  r.  In  limiting  our  interval  of  observation  to  t  instead  of  Tq  we  incur  a  slightly  greater 
ambiguity  in  the  estimation  of  namely  [see  Eq.  (51)]  , 


w' 


(63) 


We  define  6  to  be  the  output  of  our  noisy  zero-crossing  counter  at  time  t,  viz 


8  =  ?p.  n  .  (64) 

Furthermore,  since  a)'  is  bounded  between  u>M  and  u>M+1  [see  Eq.  (61)]  and  since  O  is  in  the 
same  interval  (of  width  o>M+1  -  o>M  =  27r/TQ),  we  conclude  [with  the  application  of  Eqs.  (63)  and 
(64)]  that  the  maximum  error  in  the  estimate  (9)  of  0  is 

|S-eU<„(^A)  .  <fc5) 

'  o  ’ 

Consequently,  we  would  like  to  make  r  and  Tq  as  large  as  possible;  however,  note  that  at  the 
times  of  interest  (namely  at  the  peaks),  (e(l/2,  r)  ]  behaves  essentially  as  cos(ttt/To).  Thus,  as 
r  is  increased,  the  condition  expressed  in  Eq.  (62)  becomes  more  difficult  to  satisfy,  therefore  n( 
must  increase  (see  below).  This  represents  a  trade  among  nQ,  Tq,  and  r.  To  be  specific,  we 
choose  a  reasonable  value  for  r,  namely,  r  =  Tq/3.  Then,  from  Eq.  (59),  |e(l/2,  Tq/3)  |  >  l/4 
and  we  obtain  the  following!  from  Eqs.  (62)  and  (65), 


1-cn(0) 

cn(0) 


>8 


(66) 


t  Here  we  ore  neglecting  the  restriction  that  t;  in  so  doing,  we  make  a  slight  error 
insignificant  for  lorge  M  (M^-10,  typicoliy). 

♦  Equotion  (67)  therefore  requires  5tr/!0  4  P  See  Eq.  (46) 


However,  this  error  is 


t5 


Note  that  Eq.  (66)  specifies  an  upper  bound  for  cn(0);  we  may  therefore  use  this  as  the  setting  for 
the  threshold  in  our  noisy  zero-crossing  counter.  That  is,  6  =  l/9. 

We  now  consider  the  behavior  of  cn(0).  Specifically,  cn(t)  is  the  transform  of  all  spectral 
lines  except  the  two  pairs  at  to  =  ±caM,  ±M+1-  Therefore,  cn(0)  is  merely  the  sum  of  the  mag¬ 
nitudes  of  all  such  lines,  viz., 


cn(0)  =  l  -  S<P>«V  -  S<P><  -  «M)  -  S^(o,M+  ,)  -  S<P><  -«M+1) 


(68) 


We  now 


consider  the  worst^  case,  namely. 


=  Sn 


Since  s'P»tom)  ■  S>„  ),  we  observe*  that  the  coefficient  Cfi[as  defined  by  Eq.  (66)]  is 
Cn=4SnP><"M>  ' 


OO  /  D  \ 

Furthermore,  since  S  sj^u^)  =  1  [see  Eq.  (32)]  ,  we  may  rewrite  Eq.  (68)  as 


n(0)=  1  — 4S^P)(o,M)=  1-Cn 


(69) 


The  behavior  of  S^c^)  [recall  that  sf  >(caM)  >  S^m)  for  all  m  *  *M,  *(M  +  D1  may  be 
determined  by  considering  the  extension  of  Eq.  (19)  to  periodic  functions,  vi£., 


c<P),„  ) 

n  aM 


>/pi<"m»2 


*  '4p,«“m>|2 


m=  -« 

For  convenience,  we  temporarily  adopt  the  notation. 


h  =  S{P)(w  )  . 

mom 


Thus,  Eq.  (69)  becomes 


4h 


cn(0)  =  1 


M 


(70) 


E  h 


We  adopt  the  further  notation. 


f  That  is,  the  case  which  requires  the  largest  n^  This  corresponds  to  P  -  1/2.  ^ 

♦  Note  that  the  transform  af  any  corresponding  pair  of  spectral  lines,  say  5^  (uj  ond  15  mere*y 

2S<P)(u  )  cos u  t,  since  oil  components  of  are  rea'  onc*  passive. 


Equation  (70)  then  may  be  written  as 


cn(0)  =  l- 


4h:M  +  s' 

M  m  m 


(p) 

We  recognize  that  cn(0)  (and  therefore,  nQ)  is  maximized  when,  for  fixed  hM  =  So  ("M)*  the 

sum  S'  h2"  is  maximized,  In  order  to  find  the  maximum  of  this  sum,  we  find  it  convenient  to 
m 

state  Tnd  solve  an  equivalent  problem,  as  stated  in  Lemma  2  below.  In  this  lemma  we  recognize 
that  the  quantities  Sbk<  at  and  N  correspond  respectively  to  S'hm>  1  -  4hM  and  2n.  The  value 


of  a2  will  become  apparent  shortly. 


Lemma  2.' 

Given  (bk)  k  =  1, 2,  3, 


such  that 


0  ^  a2 


Then  the  set  {b*}  which  maximizes 


V  wN 

L  \ 


k  =  1,2, _ K-  1 


—  (K  —  1)  a2  k  =  K 


and  [x]  =  the  maximum  integer  contained  in  x. 

Maximizing  c  (0)  corresponds  to  maximizing  nQ  lie..  it  makes  the  mequultts  t  1 
difficult  to  satisfy"  This  "worst"  condition  as  described  in  Lemma  2  may  be  combm.  « 
Eq.  (71)  to  give 


fThe  proof  of  Lemmo  2  will  be  found  in  the  Appendix. 
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(73) 


4h 


c(0)4  1  - 


M 


4hM+(K-1)a2  +{a!  -  <K  -  1)  a2) 


We  recognize  that 


al  _  1  4hM 


K  -  1 


1  -  4hM  | 

a2  1 


As  for  a2>  which  represents  the  maximum  value  of  (u>m)  for  m  ^  ±*M  +  1^’ 

back  to  Eq.  (53)  and  to  the  discussion  preceeding  Eq.  (54)  to  obtain 

2 


a2  "  hM 


2  ax 


A  portion  of  the  denominator  of  Eq.  (73)  may  be  bounded  above  as  follows: 


(K  -  1)  a2 


n  ?n  ?nl(1-4h\ll  /1_4hM  I1  4hMlV  | 

+  {ar(K-i,a2}2  a2  |  +  "  l“l)  I 


"  ^  1  ~  4hM  j 


a. 


Thus,  Eq.  (73)  becomes 


Cn<0)  ^  1  " 


/  1  ~  4hM\  2\Z 

\  4a 2  '  'hM' 


1  + 


Recalling  the  definition 


2air 

e0  =  rr.3 


we  substitute  Eq.  (74)  in  Eq.  (75)  to  obtain 


cn(0)«  1 


/  1  ~  4hM  W  nM  *o\ 

hM  ; 


i  + 


l4h„„  -  4e 

M  o 


h„  -  e_ \ 2 

M 


We  now  consider  that  value  of  hM  (HM>  let  us  say)  which  maximizes 


E  = 


4hu  -  4c 
M  o 


[  M  €oY 
'  hM  > 


-W=-‘)  (’-4) 


->n  A 
€  >2  -1 


Differentiating,  we  get 


refer 


(74) 


(75) 


(76) 


(77) 


(78) 
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dE 


_ _ L-  U  -  +  (2n  -  1)  - l)  (l  -  T-2-)  (~2~) 

4h^j  V  hM '  '4hM  '  '  hM '  hM 


€  vZ  -2/  6 


dh 


M 


Setting  this  derivative  equal  to  zero,  we  find  that 


2ne 


HM  = 


1  +  4eo(2"  -  1) 


We  note  that  e_  ^  HM  4:  l/4  as  of  course  it  must.  Substituting  the  value  hM  =  HM  into  Eq.  (78)  yields 


'  o 

_  1  Z1  ~  4eo\  (i  _  2-n  -  4e  +  4e  2_n)2  _1 

"  4  V  2ne  /  °  ° 


2"n  2n 

*h~  (1-4eo> 

o 


Thus,  Eq.  (77)  becomes 


cn(0)^  1  - 


2  -n  2n 

o 


(79) 


Returning  to  the  condition  expressed  in  Eq.  (66),  and  using  Eq.  (79)  we  obtain  as  our  condition, 
after  some  algebra. 


l4  .  ,2%  fo2 

(1  ~  4eo>  ^  ~T~ 


(80) 


where  e  is  defined  in  Eq.  (76).  Rewriting  Eq.  (80)  in  terms  of  logarithms,  we  obtain 
o 

n  +  2n  log  J;-  >loe  Z-f  ■  <81) 

o  o 

However,  due  to  the  convexity  of  the  function  log  t/(l  -  x)  in  the  interval  0  <  x  <  1,  we  may  bound 
log  l/(l  -  x)  from  below  by  its  tangent  at  x  =  0,  viz.. 


lo®2  1  -  x  ^  [(  dx  l0g2  1  -  x*  x=  J 


log2  T=H  loS2e  ■ 

Applying  this  inequality  to  Eq.  (81)  results  in  a  slightly  larger  value  of  nQ.  Therefore,  we  obtain 

(82) 


n  +  2n+2  e  log2e 

o 


We  now  define 
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M«>1 


(P) 

Fig.  6.  Power  spectrum  Sn(u)  end  one  period  of  the  iteroted  outocorrelotion  function  (t) 
os  a  function  of  the  number  of  iterations  n. 
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and  so 


y  =  i°g2  #- 

O 


€ 

O 


z^y 


Eq.  (82)  then  becomes 

n  +  2n+3_y  log2  e  >y 

The  minimum  value  of  n  (n  ,  for  instance)  which  satisfies  Eq.  (83)  is  clearly  nQ  =  y  -  2, 


no  =  l0g2  TZ  -  2  "  l0g2  27- 


(83) 


(84) 


Thus  n  represents  the  number  of  iterations  of  the  normalized  autocorrelation  of  f  (t)  which  we 
'  O  e  P 

need  in  order  to  obtain  a  number  0  satisfying  Eq.  (67).  This  completes  the  proof  of  Theorem  2. 


VI.  EXPERIMENTAL  RESULTS 

The  procedure  described  in  Sec.  V  for  determining  the  energy  peak  of  a  finite  duration  signal 
to  a  finite  resolution  was  simulated  on  a  digital  computer.^  Results  of  the  experimentation  are 
shown  in  Fig.  6. 

Specifically,  the  signal  f(t)  chosen  was  a  small  segment  (T  =  0.0256  sec)  of  human  speech 
sampled  at  a  10  kc  rate.  Figure  6  shows  only  one  period,  the  interval  0  ^  t  <  Tq,  of  fp(t) 
and  R<P,(t)  for  n  =  0,  1,  2,  as  well  as  S^w)  for  n  =  0,  *,  2,  3.  Tq  was  chosen  equal  to  2T  for  this 
experiment.  Note  that  Sn(w)  is  shown  only  as  a  visual  aid;  its  calculation  was  clearly  not  nec¬ 
essary  for  the  generation  of  R<P,<t).  We  chose  to  show  S^d)  rather  than  S^o^)  for  convenience 
of  programming;  as  a  result,  we  observe  the  {[sin  a  fu)  —  u>m)]/a(w  —  envelope  quite  dis¬ 

tinctly  in  S  (a>). 

One  clearly  observes  the  rapid  convergence  of  both  r|^  \t)  and  S^fw)  to  the  frequency  of 
maximum  energy  density. 


VO.  APPLICATIONS  AND  CONCLUSIONS 

Theorem  1  expresses  the  fundamental  result  that  the  limit  of  the  iterated  normalized  auto¬ 
correlation  function  of  a  signal  is  a  cosine  wave  at  a  frequency  (©)  corresponding  to  the  maximum 
energy  peak  of  the  signal’s  spectrum.  However,  two  aspects  of  the  procedure  by  which^one 
arrives  at  this  limit  require  unbounded  complexity:  first,  the  interval  over  which  the  n  iter¬ 
ation  must  be  calculated  (assuming  the  signal  to  be  of  finite  duration  T)  is  |t|  ^  2  T;  second, 
the  number  of  iterations  grows  without  bound.  These  two  difficulties  require  unlimited  equip¬ 
ment  and  time,  respectively.  Clearly,  the  reason  behind  these  infinite  operations  is  that  we 
are  asking  for  an  absolutely  perfect  determination  of  ©.  Naturally,  we  are  willing  to  accept 
some  error  in  this  determination  in  any  practical  situation.  Taking  advantage  of  this  fact,  we 
are  able  to  establish  Theorem  2  in  which  we  offer  a  procedure  for  estimating  S  (to  within  an 
arbitrarily  small,  but  finite,  error)  which  requires  a  finite  number  of  calculations  over  a  fixed 
time  interval.  Thus,  by  accepting  an  error  in  the  determination  of  0,  we  have  been  able  to 
eliminate  both  undesirable  aspects  of  the  original  procedure. 

t  A  program  for  simulating  the  iterated  autocorrelation  written  by  the  author,  and  a  spectrum  analyses  program 
written  by  C.  Rader,  were  run  on  Lincoln  Laboratory’s  TX-2  high-speed  digital  computer  (see  Ref.  3). 


2i 


The  procedure  for  obtaining  the  estimate  (@)  of  0  may  be  mechanized  as  follows  (see  F'ig.  7). 
The  signal  f(t)  would  be  stored  on  a  tapped  delay  line  (of  Tq  seconds).  At  time  Tq,  the  impulse 
response  of  a  linear  filter1-  would  be  set  equal  to  the  values  of  the  taps  and  the  output  of  the  delay 
line  would  then  feed  into  the  filter.  The  output  of  the  filter  during  the  interval  2Tq  <  t  <  3Tq 
(now  equal  to  the  autocorrelation  of  the  signal)  would  then  be  fed  back  into  the  delay  line.  The 
procedure  is  repeated  nQ  times,  after  which  the  output  of  the  filter  is  sent  through  the  noisy 
aero-crossing  counter,  which  provides  the  estimate  e. 


Fi«.  7.  tmplementotion  ot  the  procedure  tor  detecting  the  frequency  of  maximum  energy  density. 


The  applications  of  these  results  appear  to  be  numerous.  The  detection  of  the  energy  peak 
of  an  arbitrary  signal  (for  example,  the  energy  peak  corresponding  to  a  time  series)  is  often  of 
great  interest,  and  may  be  detected  as  above.  Furthermore,  this  method  may  be  used  for  se¬ 
lecting  that  one  out  of  N  possible  signals  transmitted  over  a  communication  link;  in  this  appli¬ 
cation  it  is  necessary,  of  course,  that  the  location  of  the  absolute  peaks  in  the  spectra  of  the 
N  signals  be  distinct,  one  from  the  other,  such  as  is  the  case  in  Frequency  Shift  Keying.  Another 
application  may  be  found  in  locating  the  peak  frequency  in  the  return  signal  from  a  Doppler  radar. 
These  are  but  a  few  of  the  possible  applications  of  the  procedure  described. 

The  main  conclusion  to  be  drawn  from  this  Study  is  that  the  iterated  autocorrelation  pro¬ 
cedure  represents  a  new  method  for  detecting  the  frequency  of  maximum  energy  density  of  a 
signal.  Some  areas  of  application  have  been  suggested  above  briefly;  but,  more  careful  con¬ 
sideration  and  analysis  must  be  undertaken  before  one  can  determine  the  advantages  or  disad¬ 
vantages  of  this  system  compared  to  any  other.  In  fact,  one  of  the  main  purposes  of  presenting 
this  material  is  to  stimulate  thinking  about  the  possible  applications  and  merits  of  this  new  de¬ 
tection  approach. 


t  This  fitter  is  2T.  sec  tong,  with  it*  impulse  response  h(t)  adjusted  such  that  h(t  +  T0)  -  h(t),  thus  rePr«en 

the  periodic  version  of  the  signal.  In  Fig.  7,  this  linear  filter  is  represented  as  o  combination  of  *  holding ,  circuit 

and  a  second  tapped  delay  line,  the  outputs  of  which  are  multiplied  and  summed  to  form  the  output  of  the 


filter. 
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APPENDIX 

PROOFS  OF  LEMMAS 


Lemma  1. 

We  may  express  e(/3,  t)  as  the  real  part  of  the  sum  of  two  vectors  in  the  complex  plane,  viz. , 
e(/3,t)  =  Re  {/3  exp  [j(2jrt/To>  M]  +  (1  -0)  exp  [j(2Tt/TQ)  (M  +  1)]} 

=  Re  [A(/3,t)  ejY<^'t)1 

where  A(0,t)  is  the  magnitude  of  this  vector  sum,  and  y(0,t)  is  its  phase  angle.  We  may  express 
A(0,t)  and  y(0,t),  after  some  trigonometric  simplifications,  as 

,  2  2  >Ttll/2 

A(0,t)  ■=  0  +  (1  ~/3)  +  20(1  -0)  cos  ^  I  (A- 1) 

1  o 

and 


( 1  —  (3 )  sin  —■ 

y(0 ,  t)  =  M  +  tan"  1 - ~-r- 

o  0  +  (1  -0)  cos^ 


(A-2) 


Forming  dA(0,t)/d0  =  0,  we  easily  show  that 
A(j,t)<  A(0,t) 

Furthermore,  we  observe  that  for  0  <  t  ^  T  /2, 


(A-3) 


<  0 
at 


from  which  we  conclude  that 


A(0,t1)  >  A(fl,t2)  for  tt  t2«  — 


(A-4) 


Let  uk  be  that  instant  when  y(0,uk)  =  (k  —  1)  ir.  Now,  for  0 t (Tq/2),  we  observe  that  dy(/3,  t)/dt 
is  always  positive;  furthermore,  since  A(0,t)  is  a  decreasing  function  of  t  in  this  interval  |see 
Eq.  (A-4)],  we  note  that  de(0,t),/dt  passes  through  zero  each  time  the  vector  A(0,t)  e  '  passes 
through  either  the  second  or  fourth  quadrants  of  the  complex  plane.  We  may  then  conclude  that 


the  instants  u^  and  must  alternate,  viz.. 


<  uk_i  ^  Sk.<  uk  <  Sk+1  $ 


In  addition. 


I  e<0 ,  Sk)  |  ^  |  e<0 ,  uk)  |  =  A(0 ,  uk> 


A-5) 


which  follows  from  the  definition  of  and  u^  and  the  fact  that  the  relative  maxima  and  minima 
of  e(0 ,  t)  occur  when  the  vector  A<0 ,  t)  ^  passes  through  the  second  and  fourth  quadrants. 

It  is  also  clear  that 

|  e(0,  t)  |  =  |  A<0 ,  t)  cosy(0,  t)l  «  A(0,  t) 


Z3 


From  Eqs.  (A-3)  through  (A-6)  we  may  write 


|eO,S.)|  »|e(0,u  )|  =AO.uk)>A(|-,uk)>A(|,T)^|A(|-T)cosy(|,r)|  =  |e(y,r)|  , 


where 


°^Sk<uk^T^—  ■ 

This  completes  the  proof  of  Lemma  1. 

Lemma  2. 

Consider  a  sequence  of  sets  (b^)where  we  define 


fb'0'}  ■  (bk) 


Corresponding  to  each  such  sequence,  we  define 

oo 

V  rh(i)iN 
ffi  -  L  ibk  1 

k=  1 

where,  obviously,  %  =  a.  Since  we  have  freedom  in  labelling  the  subscript  k,  we  choose  to 
arrange  the  bk0)  in  a  nonincreasing  sequence;  that  is, 

b(0)  >  b(0) 
bk  ?  k+1 

The  iteration  on  i  may  be  described  essentially  as  follows.  We  concentrate  on  two  special  values 
of  k,  say  k^i)  and  k2(i),  where 

k^(i)  =  the  smallest  value  of  k  for  which  bk  <  , 

k  (i)  =  the  smallest  value  of  k  greater  than  or  equal  to  K,  for  which 

2  >  0  • 
k 

Corresponding  to  these  values,  we  define 

Di  =  bk^(i) 
di  =  bk2(i)  ' 

Then,  the  set  {b^i+1)}  is  defined  in  terms  of  {b£l}  }  by 

f  D.  +  a  .d.  for  k  =  k.(i) 

I  l  ii  1 


bii+1)  ={di(1-“i 


(1  —  as)  for  k  =  k-,(i) 


otherwise 
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/  a,  —  D  .  v 

'i  =  ^  l1'  - d~) 


Thus,  at  each  stage,  we  transfer  a  fraction  (a.)  of  to  being  careful  not  to  allow  any 

to  exceed  a2-  Now,  it  is  clear,  from  the  definition  of  K  and  Eq.  (A-7),  that  k^i)  <  K  and 

k  (i)  >  K  For  those  i  where  k,(i)  remains  constant,  D  is  an  increasing  function,  and  for  those 
2  1  1  im  im 


i  where  k2(i)  remains  constant,  d^  is  a  decreasing  function.  Also,  since 


b,(0)  ^  b}°\  we  have 


D.  >.  d.  for  all  i 
i  l 


Furthermore,  at  each  step,  we  conserve  the  sum 


l  ^  =  *1  • 

k=  1 

The  iterative  procedure  continues  until  we  reach  that  value  of  i  (i^,  lt?t  us  say)  for  which 

<i> 

b,  °  =  a,  k  1,  2,  ...  ,K  -  1  . 

k  2 

We  now  show  that  +  ^  ^  for  all  i  <  iQ.  As  a  result  of  Eq.  (A-7)  we  may  write 

<Vi  =  "i  -  °iN  -  diN  +  (Di +  aidi)N  + 1(1  -  ai>  di]N  • 

Using  the  binomial  expansion,  we  get 


2  (")  iDi'“idi|N'J  *  dii<-«id!iN'ii-DiN-di 

J  =  0 


=  Z  (y)(aidi)N’j  [Di +  (-1)N_J  di )  • 

3  =  0 


But,  since  ^  d.,  we  have,  for  i  <  iQ, 


°-i+l  -  erf  >  ° 


We  now  consider  the  remainder 


l  \°  =a1-(K-l)a2HV  . 


i  Cf 

k  =  K  <•  3  lk=K  1 


and  this  upper  bound  may  be  achieved  when 
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‘V 


y  k  =  K 
0  k  >  K 


We  observe  that 


{\°] )  =  ’ 


(A-9) 


where  {b,|lo>}  is  defined  by  Eqs.  (A-8)  and  (A-9).  Thus,  we  have  shown  that  starting  from  any 
{b^0*},  the  sequence  {b^1*}  may  be  defined  in  a  way  such  that 

lim  {b^1*}  =  {b*} 

and 

CTi <  CTi+l  i<:io 

This  completes  the  proof  of  Lemma  2. 
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